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ABSTRACT 

We report evidence that the gas properties of X-ray groups containing radio galaxies differ 
from those of radio-quiet groups. For a well-studied sample of 7?<9&4r-observed groups, we 
found that more than half of the elliptical-dominated groups can be considered "radio-loud", 
and that radio-loud groups are likely to be hotter at a given X-ray luminosity than radio-quiet 
groups. We tested three different models for the origin of the effect and conclude that radio- 
source heating is the most likely explanation. We found several examples of groups where 
there is strong evidence from Chandra or XMM-Newton images for interactions between the 
radio source and the group gas. A variety of radio-source heating processes are important, 
including shock-heating by young sources and gentler heating by larger sources. The heating 
effects can be longer-lasting than the radio emission. We show that the sample of X-ray groups 
used in our study is not significantly biased in the fraction of radio-loud groups that it contains. 
This allows us to conclude that the energy per particle that low -power radio galaxies can inject 
over the group lifetime is comparable to the requirements of structure formation models. 
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1 INTRODUCTION 

Radio galaxies must be transferring large quantities of energy to 
the surrounding group- or cluster-scale gas. The PdV work done 
on the gas by source expansi on is of the order of 10 52 - 10 53 J for 
source ages ~ 10 8 years (e.g. lCroston et all2003l) . The first strong 
evidence for radio-source heating was recently found in the nearest 
radio galaxy, Centaurus A (Kraftetal. 2003): Chandra observa- 
tions reveal a prominent shell of emission capping the inner south- 
western radio lobe, which has a temperature ten times that of the 
surrounding galactic atmosphere, providing strong evidence that 
Cen A is shock-heating its atmosphere. There is also evide nce for 
heatin g in the powerful FRII radio galaxy, Cygnus A l Smi th et all 
120021) . where a slight temperature increase in the X-ray gas is seen 
at the front edges of the lobe cavities. Another recent result helps to 
emphasise that more than one energy transfer mechanism is likely 
to operate. Deep Chandra observations I F abian et all20 03j) of the 
Perse us cluster [t he first cluster where "cavities" were observed 
( Bohri nger et alj Fl 993)1 have revealed the presence of ripples in 
the cluster gas that are suggestive of sound waves emanating from 
the central radio source, 3C 84. Fabian et al. calculate that the time 
scale between successive wavefronts is comparable to estimates of 
the radio-source lifetime, so that it seems plausible that an intermit- 
tent radio source is producing the observed ripples. 

There is evidence that energy injection of some sort is required 
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to explain the observational properties of X-ray groups and clus- 
ters, which do not agree with the prediction s of CDM models of 
struct ure formation (e.g. lArnaud and Evrardlll999l : IPonman et alJ 
1999). Several authors have considered th e effects of heatin 
due to feedback from star formation (e.g. Balog h"et al J 1 199' 
iBrighentian d Mathews 2001). However, in most cases it is found 
that heating from star formation and supernovae does not provide 
sufficien t energy to exp lain the prop erties of both gro ups and clus- 
ters (e.g. lKravtsov and YepeJ200Cl) . IWu et all <2000t) find that su- 
pernova heating can generate only ~ 1/10 of the energy required. 
These results have le d m any authors to c onsider AGN heating. 
IChurazov et alJ 1200 ll) and lBohringer et all 120021) consider the ef- 
fects of radio-galaxy heating on cluster structure and conclude that 
sufficient heating could be provided. Evidence from recent ob- 
servations of clusters, as well as simulations, suggest that AGN 
could provide a sufficie ntly distributed heating mechanism for 
this m ethod to work (e.g. iBriiggen and KaisejEo02t iFabian et alJ 
2003), which removes a long-standing objection to radio-source 
heating models. Since radiative cooling in the absence of feedback 
overpredicts the mass in galaxies (e.g. Col a 199 11) . it seems likely 
that some sort of feedback is required to explain observed c luster 
properties and the galaxy luminosity function (e.g. lVoit and Brvar] 
1200 ltlBenson et all2003l : lKavl2004l) . 

In addition, nearly all cooling-flow clusters contain a central 
radio galaxy (e.g. lEilekll2004l) . and there is considerable evidence 
that the radio galaxy displaces gas in the cooling-flow regions (e.g. 
iBohringer et alll993l : lBlanton et alfeOollFabian et all2003l) . It is 
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therefore important to consider the different means by which ra- 
dio galaxies could influence the observed properties of the cool- 
ing flow, and in particular whether they can help explain why large 
quantities of gas do not ap pear to cool past ~ 1/3 the outer c luster 
temperature (e.g. IPeterson et al J 120031 ISakelliou et aljEo02l) . Al- 
though non-heating solutions to this problem exist, it is most plau- 
sible that the bulk of the gas in these systems is reheated by one 
of several possible mechani sms, such as cluster mergers , thermal 
conducti on (e.g.lVoigt et all2002t|Voigt and Fabiarl20o4). o r AGN 
heating jiinnev and TaboJll995tlBruggen and Kaisej200lL l2002t 
iRevnolds et al 120 02): the last of these is particularly attractive be- 
cause of the possibility of self-regulation via feedback. Investigat- 
ing how radio galaxies can affect their hot-gas environments is 
therefore important to our understanding of several problems of 
cluster physics. 

In lCroston et alji2003h . we showed that radio galaxies have an 
important impact on groups, and found evidence that the properties 
of "radio-loud" and "radio-quiet" groups differ, in the sense that 
"radio-loud" groups are hotter than "radio-quiet" groups of com- 
parable X-ray luminosity. That work used a fairly small and inho- 
mogeneous sample, and the analysis methods were comparatively 
basic. Here we present a detailed analysis of a larger, homogeneous 
sample of groups (Osmond and Ponman 2004) in order to confirm 
and investigate further the conclusions of the earlier work. 



2 SAMPLE SELECTION AND ANALYSIS 
2.1 The elliptical-dominated sample 

This analysis uses the GEMS group sampl e who se X-ray properties 
were presented bv lOsmond and Ponmanl J2004) (hereafter OP04). 
The sample consists of 60 groups, including loose and compact 
groups that may be spiral- or elliptical-dominated. For the larger 
GEMS sample, the LxlTx relation shows mo re scatter than for 
the earlier work of Helsdon and Ponman (2000), so that it may be 
more difficult to distinguish heating effects. 

We considered only the subset of groups in the OP04 sample 
that contain a large elliptical galaxy, as spiral-dominated groups are 
unlikely to possess a strong radio galaxy. We found that the scatter 
in the LxjTx relation is significantly reduced in our subsample, 
suggesting that spiral-dominated groups have different gas prop- 
erties from those with a large elliptical galaxy. OP04 did not re- 
port a significant difference in the Lx /Tx properties of elliptical- 
and spiral-dominated groups; however, they do find that most X-ray 
bright groups contain a bright central early-type galaxy. They used 
the group's spiral fraction to compare group properties. Our classi- 
fication on the basis of the dominant galaxy's morphology may be a 
more useful measure of the group's history and current properties. 
This is supported by OP04's conclusions about the importance of 
elliptical brightest group galaxies. 

In order to determine whether each group contains a large el- 
liptical galaxy, we downloaded DSS2 (Digitized Sky Survey 1 ) im- 
ages of each group. Typically groups were either dominated by one 
large galaxy, or else there were several bright galaxies of similar 
magnitude. In the first case, we rejected any groups whose dom- 
inant galaxy is a spiral or SO. In the second case, where there 
was no obviously dominant galaxy, the group was only rejected 
if none of the bright galaxies is an elliptical. In some cases it 
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was not possible to tell by eye whether a galaxy has an ellipti- 
cal or SO morphology, and so we followed up all of the elliptical 
groups using Simbad and NED to confirm the galaxy morphol- 
ogy. Unfortunately there are many cases where these databases dis- 
agree, and multiple classifications exist in the literature. We there- 
fore only rejected groups where the ambiguous galaxy was clas- 
sified as SO by both Simbad and NED. In all we excluded 15 
spiral-dominated groups: HCG 10, HCG 15, HCG 16, HCG 40, 
HCG 68, HCG 92, NGC 1332, NGC 2563, NGC 3227, NGC 3396, 
NGC 4565, NGC 4725, NGC 5689, NGC 5907 and NGC 6574. In 
the course of following up galaxy morphologies, we also found one 
group, HCG 22, for which the supposed group members covered an 
implausibly large range in redshift. As it is unclear which galaxies 
are real members of this group, we excluded it as well. 

In addition to excluding spiral-dominated groups, we also 
had to exclude 10 groups for which OP04 could not make X-ray 
spectral measurements: HCG 4, HCG 48, HCG 58, NGC 1808, 
NGC 3640, NGC 3783, NGC 4151, NGC 4193, NGC 6338 and 
NGC 7714. As the aim of the study was to compare the properties 
of radio-loud and radio-quiet groups, we also had to exclude two 
groups with low declinations that are not covered by the surveys 
used to identify the radio sources, NGC 1566 and NGC 7144 . Fi- 
nally, NGC 315 was excluded, because Iworrall et alj J2003T) find 
a dominant contribution to the X-ray luminosity from the AGN 
and X-ray jet, and HCG 67 was excluded because of ambiguity 
in whether or not an identified nearby radio source is actually as- 
sociated with the group. The final sample of elliptical-dominated 
groups contained 30 members. In TableQ we list the groups in the 
sample along with their redshift and the X-ray properties measured 
by OP04 used in this work. 

2.2 Identifying radio sources 

We then divided the elliptical-dominated groups into radio-loud 
and radio-quiet subsamples based on the properties of any ra- 
dio source associated with eac h group. For each group in the 
sample, we used N ED, NVSS ICondon et all 1 19981) and FIRST 
iBecker et ai]ll995h to locate radio sources that could potentially 
be associated with a group member. We first checked NED for any 
known radio galaxies associated with the group's galaxies. If an as- 
sociated radio source was found, we adopted the 1.4-GHz radio flux 
from the NED database. For any groups where a radio source was 
not found in this way, we searched for sources within a radius of 10 
arcmin using NVSS and FIRST. We then checked the location of 
each candidate for an associated radio source on the DSS2 optical 
images to ensure that the radio source was roughly at the centre of 
a group galaxy (using SIMBAD to confirm that the galaxy is at the 
redshift of the group). If this was not the case, then we rejected the 
radio source on the basis that it was probably a background object. 
This method should be a more robust way of defining radio-loud 
and ra dio-quiet groups than the method we used in lCroston et alJ 
(2003), as the latter may contain spurious "radio-loud" designa- 
tions if any of the radio sources were not in fact associated with the 
group. Table [2] gives the 1.4-GHz radio flux and luminosity den- 
sities, Li.4, and the location for all of the associated sources. We 
also list the distance between the radio source and group centre, de- 
fined in OP04 as the position of the group-member galaxy nearest 
to the centroid of the X-ray emission. In two cases (NGC 5171 and 
HCG 90) the radio source is a significant distance from the group 
centre, which may reduce the likelihood that it could strongly affect 
the group gas properties; however, as in both cases the radio-source 
is clearly associated with a group galaxy and lies in the extended X- 



Evidence for radio-source heating of groups 3 

Table 1. The elliptical-dominated groups sample. Properties listed here are taken from OP04. /3-model parameters marked with a star are estimated using the 
method described in the text. 



Group 


Redshift 


Temperature (keV) 


Abundance (Zq) 


log(L x ergs 1 ) 


Tcut (kpc) 




r c (kpc) 


HCG 42 


0.0128 


0.75±0.04 


0.29±0.10 


41.99±0.02 


112 


0.56 


4.69 


HCG 62 


0.0146 


1.43±0.08 


2.00±0.56 


43.14±0.04 


282 


0.48 


2.44 


HCG 90 


0.00880 


0.46±0.06 


0.08±0.03 


41.49±0.05 


101 


0.41 


0.91 


HCG 97 


0.0218 


0.82±0.06 


0.23±0.10 


42.37±0.05 


339 


0.44 


2.73 


NGC 383 


0.0173 


1.51 ±0.06 


0.42±0.08 


43.07±0.01 


633 


0.36 


2.11 


NGC 524 


0.00793 


0.65±0.07 


0.22±0.15 


41.05±0.05 


56 


0.45* 


0.37* 


NGC 533 


0.0181 


1.08±0.05 


0.68±0.23 


42.67±0.03 


372 


0.42 


2.21 


NGC 720 


0.00587 


0.52±0.03 


0.18±0.02 


41.20±0.02 


65 


0.47 


1.15 


NGC 741 


0.0179 


1.21 ±0.09 


2.00±0.67 


42.44±0.06 


386 


0.44 


2.30 


NGC 1052 


0.00492 


0.41±0.15 


0.00±0.02 


40.08±0.15 


25 


0.45* 


0.04* 


NGC 1407 


0.00565 


1.02±0.04 


0.23±0.05 


41.69±0.02 


105 


0.46 


0.08 


NGC 1587 


0.0122 


0.96±0.17 


0.47±1.24 


41.18±0.09 


77 


0.46 


4.34 


NGC 3557 


0.00878 


0.24±0.02 


0.00±0.01 


42.04±0.04 


95 


0.52 


1.13 


NGC 3665 


0.00692 


0.47 ±0.10 


0.17±0.14 


41.11±0.08 


71 


0.47 


1.08 


NGC 3607 


0.00411 


0.35±0.04 


0.23±0.10 


41.05±0.05 


62 


0.39 


1.98 


NGC 3923 


0.00459 


0.52±0.03 


0.18±0.05 


40.98±0.02 


34 


0.55 


0.63 


NGC 4065 


0.0235 


1.22±0.08 


0.97±0.48 


42.64±0.05 


425 


0.36 


3.08 


NGC 4073 


0.0204 


1.52±0.09 


0.70±0.15 


43.41±0.02 


470 


0.43 


9.42 


NGC 4261 


0.00785 


1.30±0.07 


1.23±0.42 


41.92±0.03 


112 


0.44 


40.08 


NGC 4636 


0.00565 


0.84±0.02 


0.41±0.05 


41.49±0.02 


68 


0.47 


0.30 


NGC 4325 


0.0252 


0.82±0.02 


0.50±0.08 


43.15±0.01 


307 


0.58 


27.56 


NGC 4589 


0.00676 


0.60±0.07 


0.08±0.03 


41.61±0.05 


122 


0.52 


9.33 


NGC 4697 


0.00454 


0.32±0.03 


0.07±0.02 


41.01±0.02 


53 


0.46 


1.25 


NGC 5044 


0.00820 


1.21 ±0.02 


0.69±0.06 


43.01±0.01 


180 


0.51 


5.96 


NGC 5129 


0.0232 


0.84±0.06 


0.66±0.28 


42.33±0.04 


151 


0.43 


3.14 


NGC 5171 


0.0232 


1.07±0.09 


1.47±1.25 


42.38±0.06 


298 


0.45* 


81.26 


NGC 5322 


0.00702 


0.23±0.07 


0.00±0.02 


40.71±0.10 


43 


0.45* 


0.18 


NGC 5846 


0.0063 


0.73±0.02 


1.25±0.69 


41.90±0.02 


94 


0.51 


2.19 


NGC 5930 


0.00969 


0.97±0.27 


0.17±0.12 


40.73±0.07 


29 


0.45* 


0.19* 


IC 1459 


0.00569 


0.39±0.04 


0.04±0.01 


41.28±0.04 


121 


0.45 


0.74 



Table 2. Radio sources associated with groups. Coordinates are for the radio-source position. Column 3 (D) is the physical distance between the radio source 
and the group centre (see text for discussion). Luminosity densities are determined by assuming that the source is at the redshift of the group. 



Group 


RA (J2000) 


Dec (J2000) 


D/kpc 


Fi. 4 /Jy 


L1.4AVHZ- 1 


Ref 


NGC 383 


01 07 24.9 


+32 24 45 


5.0 





3.4 xlO 24 


NED 


NGC 524 


01 24 47.7 


+09 32 21.7 


0.51 


0.0035 


5.4 xlO 20 


NVSS 


NGC 533 


01 25 31.3 


+01 45 33 





0.0291 


2.2 xlO 22 


NED 


NGC 741 


01 56 22.1 


+05 37 39.8 


6.2 


1.0 


7.2 xlO 23 


NVSS 


NGC 1052 


02 41 04.8 


-08 15 21.1 


0.01 


0.913 


5.1 xlO 22 


NVSS 


NGC 1407 


03 40 11.9 


-18 34 49.0 


0.20 


0.0877 


6.8 xlO 21 


NVSS 


NGC 1587 


04 30 39.9 


+00 39 42.1 


0.22 


0.123 


4.1 xlO 22 


NVSS 


NGC 3557 


11 09 57.7 


-37 32 20.7 


0.91 


0.484 


1.6 xlO 23 


NVSS 


NGC 3607 


11 17 01.6 


+18 08 44.8 


30.0 


0.0181 


6.4 x 10 20 


NVSS 


NGC 3665 


11 24 43.4 


+38 45 44 





0.113 


1.2 xlO 22 


NED 


NGC 3923 


11 51 05.0 


-28 46 10.5 


2.3 


0.0312 


1.7xl0 21 


NVSS 


NGC 4261 


12 19 23.2 


+05 49 31 





18.0 


2.0 x 10 24 


NED 


NGC 4636 


12 42 50.3 


+02 41 18.7 


0.64 


0.299 


1.3 xlO 22 


NVSS 


NGC 5044 


13 15 24.0 


-16 23 7.6 


0.27 


0.0347 


5.0 xlO 21 


NVSS 


NGC 5171 


13 29 47.9 


+11 42 33.8 


186 


0.025 


3.1 xlO 22 


NVSS 


NGC 5930 


15 26 06.7 


+41 40 21.0 


3.7 


0.108 


1.9 xlO 22 


NVSS 


HCG 62 


12 53 05.6 


-09 12 21 


1.7 


0.0049 


1.6 xlO 21 


NVSS 


HCG 90 


22 02 02.0 


-31 52 10.5 


81 


0.0368 


1.4 xlO 23 


NVSS 


IC 1459 


22 57 10.7 


-36 27 43.0 


0.19 


1.28 


9.1 xlO 22 


NVSS 
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Figure 1. Plot of log(Xx) vs. log(r c ), illustrating the trend between these 
quantities. The OLS bisector fit to this line is overplotted and has a slope 
of 1.11 and an intercept of 41.53 at log 10(r c ) = 0. We used this relation 
to estimate r c for the groups where OP04 were unable to measure /3-model 
parameters. 



ray halo detected with ROSAT, we include these two radio sources 

in the sample. 

In our earlier analysis <Crostonetail2 003). we chose a single 
cut-off to discriminate between radio-quiet and radio-loud groups. 
However, at lower luminosities it is difficult to distinguish be- 
tween AGN-related radio emission and emission from other pro- 
cesses, such as starbursts. The radio source population is domi- 
nated by starbursts at 1.4-GHz luminosities below ~ 10 23 W Hz -1 
I Sadler et al. 2002). All of the radio source identifications are with 
elliptical galaxies, so that a starburst origin is unlikely, and the ma- 
jority of the sources have also been previously identified as AGN. 
However, we may detect weak AGN emission in elliptical galaxies 
that have never had a radio source sufficiently powerful to have af- 
fected the group properties. We therefore tested the importance of 
the radio-luminosity cut-off, by carrying out all of the analysis that 
follows for three choices. We first used a cut-off (cO) of Lf^ = 0, 
so that the possession of any radio source above the NVSS flux den- 
sity limit meant that a group was considered to be radio-loud. We 
then chose two higher cut-offs, based on the luminosity density of 
NGC 3665, a comparatively weak d ouble-lobed radio galaxy, as in 
the analysis of lCroston et alji2003h . These cut-offs (cl and c2) are 
L\ u l = 1.2 xlO 21 W Hz" 1 (Livccsees/lO) and Lf 4 ' = 6 xlO 21 
W Hz" 1 (LATGC3665/2). Tabled gives the total number of groups 
in each subsample for the three choices of radio-luminosity cut-off. 
We note that the NVSS flux limit of 2.3 mJy may introduce some 
bias into the selections, as this corresponds to a limiting luminosity 
of 3 x 10 21 W Hz -1 for the highest redshift group, NGC 4325, at 
z — 0.0252. This limit is close to the cut-off luminosities, so that 
for cut-offs cO and cl a few high-redshift groups could have been 
incorrectly classed as radio-quiet despite possessing a radio-source 
more luminous than the cut-off luminosity. There are five groups 
with sufficiently high redshift that a radio source more luminous 
than cl could have been missed. However, a radio source of lumi- 
nosity > c2 would be detectable in all of the groups, so that the 
results for this cut-off should provide a check for whether this bias 
is important. 



2.3 Radio sources in the OP04 parent population groups 

We found a surprisingly high fraction of radio sources in the sam- 
ple of elliptical-dominated groups from the OP04 catalogue (19/30 
= 63 per cent, assuming cut-off cO). It is often stated that radio 
galaxies are not common; a sma ll fraction of elliptical galaxies 
[e.g. ~ 5 per cent, Schmidt 1 1978)] host a large radio galaxy. How- 
ever, the fraction is cert ainly higher for the br ightest ellipticals (e.g. 
Birki nshaw and Daviesll985l) . and lHol ll99Sl) discusses the recent 
detections of small radio cores in many nearby elliptical galaxies, 
concluding that these are likely to be low-luminosity AGN. Since 
the preferred environment of radio galaxies may be the cent res of 
elliptical-dominated groups or poor clusters (e.g. lBestH2004l) . the 
high fraction of "radio-loud" groups in the sample may not be un- 
expected. If radio galaxies are found in such a high fraction of 
this type of group, then our results will have important implica- 
tions for the properties and evolution of groups. For this reason, it 
is crucial to test whether or not the OP04 X-ray observed sample 
of groups may be biased in its radio properties. The OP04 sample 
was chosen by merging nine catalogues of optical groups and then 
cross-correlating the resulting list with the ROSAT observing log. 
The parent catalogue is unbiased with respect to the groups' radio 
properties; however, it is possible that the ROSAT archive contains 
a high fraction of groups with active galaxies, since ROSAT ob- 
served many radio galaxies. This could bias the OP04 sample to- 
wards groups containing radio galaxies, although groups with pre- 
viously known radio galaxies make up a fairly small fraction of the 
OP04 sample. 

To test whether the fraction of radio-loud groups in our sam- 
ple is biased, we looked at the parent catalogues used by OP04. As 
we were only interested in the properties of elliptical-dominated 
groups, we wanted to use an electronically available catalogue 
that contained information about the morphology of the dominant 
galaxy in each group. The whole-sky group catalogue of iGarcij 
1 1993), taken from the Lyon-Meudon Extragalactic Database fits 
these criteria. It contains 485 groups having z < 0.02. This sam- 
ple is large enough that we can test whether its radio properties are 
consistent with those of the OP04 sample. 

Using Vizier 2 , we extracted from the Garcia catalogue all 
groups whose dominant galaxy has type E (elliptical) or L (SO) 
(these classifications are taken from lde Vaucouleurs et alj|l99lh . 
Although we excluded groups with a dominant galaxy with a con- 
vincing SO designation in our sample definition, we decided to in- 
clude them here, for two reasons. Firstly, as mentioned earlier, we 
found several dominant galaxies in our sample with SO designa- 
tions where later work revealed them to be misclassified ellipticals. 
Secondly, a surprisingly low fraction of the groups in the sample 
had "E" designations, so that the test sample would have been quite 
small. Including all the SO groups means that the resulting radio- 
loud fraction will be a conservative lower limit, as many of the SO 
identifications will be correct. The final sample of E and SO groups 
from the Garcia sample contains 135 groups (~ 30 per cent of the 
original sample). 

We then cross-correlated the Garcia E and SO groups with 
NVSS, searching for radio sources within 15 arcsec of the centre of 
the dominant galaxy. This method is not as accurate as the method 
we used for the OP04 radio identifications; however, the more de- 
tailed method is too cumbersome for this larger sample. To ensure 
a fair comparison, we carried out the same cross-correlation for our 
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Table 3. Samples sizes for different choices of radio-luminosity cut-off, Lf 1 ^' as given in the text. 



Cut-off number L1.4/WHZ 1 Number in RQ sample Number in RL sample 



cO 11 19 

cl 1.2 xlO 21 13 17 

c2 6xl0 21 16 14 



elliptical-dominated groups sample, using the coordinates given in 
OP04. 

For the E/SO Garcia subsample, we found radio sources asso- 
ciated with 41 groups. 32 groups had coordinates outside the re- 
gion covered by NVSS, so that the final "radio-loud" fraction of 
the Garcia subsample is 41/103, or 39.8 per cent. For our elliptical- 
dominated sample (Table 0, using the same method, we found 
a "radio-loud" fraction of 16/30, or 53.3 percent. Therefore the 
"radio-loud" fraction of the elliptical-dominated groups in OP04 
is consistent with that in this parent catalogue. The fraction we ob- 
tained here for our OP04 subsample is slightly lower than that ob- 
tained using the more accurate identification method described in 
Section l2~2l which was 19/30, or 63 percent. This is unsuprising, as 
the detailed method would find sources associated with any large 
elliptical in the group, whereas this cross-correlation method will 
only find sources associated with the dominant galaxy. 

We conclude that the OP04 group sample is not excessively bi- 
ased with respect to the groups' radio properties. The true fraction 
of elliptical-dominated groups with radio sources may be ~ 40—50 
percent (since our result for the Garcia sample is a conservative 
lower limit), rather than being as high as is suggested by our orig- 
inal analysis of the OP04 sample. It is interesting that such a high 
fraction of elliptical-dominated groups is likely to possess an AGN- 
related radio source. If all elliptical-dominated groups are capable 
of hosting radio sources, this result could indicate that radio galax- 
ies have a high duty cycle. Radio sources with obvious double- 
lobed structure make up roughly half of the radio sources in the 
sample, so that perhaps only half of the 40 - 50 percent of elliptical- 
dominated groups with radio sources could be considered to be in a 
very active state. Nonetheless, this suggests that a duty cycle where 
every elliptical-dominated group contains a radio source that is ac- 
tive for > 1/4 of the time. Weaker sources with no detected double- 
lobed structure may be in a less active stage, and the groups with no 
radio source above the NVSS limit may be in the least active phase. 
If this model is correct, then the effects of radio sources are likely to 
be important at some level in all elliptical-dominated groups. The 
analysis of the Garcia sample found that E/SO-dominated groups 
made up ~ 30 percent of the sample. Although this is a minor- 
ity of groups, E/SO-dominated group s are more likely to host a 
group-scale X-ray atmosphere lOsmond and Ponmanll2004T) . They 
are therefore likely to be in a more relaxed state, so that their X-ray 
properties will be of more relevance to structure-formation models. 



2.4 Issues with X-ray luminosity comparisons 

The analysis we carried out in ICroston et al ] <2003l) is based on 
the Lx /Tx relations for radio-loud and radio-quiet groups. In that 
analysis we did not take into account the choice of radius to which 
the X-ray luminosity was measured. OP04 used an X-ray extraction 
radius defined by the extent of X-ray emission at a significant level 
above the background. They then used the luminosities obtained 



for these regions and their fitted /3-model parameters to extrapolate 
the luminosity to a fixed overdensity radius, rsoo (= the radius cor- 
responding to 500 x the critical density of the Universe). We were 
concerned that the choice of radius used for the luminosity deter- 
mination might affect the results. There are problems in this con- 
text with both radii used by OP04. A cut-off radius defined by the 
extent of X-ray emission biases the X-ray luminosity-temperature 
(Lx/Tx) relation in the sense that a smaller fraction of the at- 
mosphere will be measured for fainter groups, because the surface 
brightness drops below the background at a smaller radius. How- 
ever, choosing a cut-off at rsoo may not be suitable for our analysis 
either, as the method used by OP04 to define rsoo is temperature- 
dependent, so that the luminosity is measured to a larger physical 
radius for hotter groups. If groups have been heated by the presence 
of a radio source, the choice of rsoo will have the effect of reducing 
the significance of any heating effect we measure. 

Clearly, it would be preferable to use the luminosity integrated 
out to infinity. Unfortunately this is not possible using a /3-model 
representation of the groups, as the solid angle integral of the sur- 
face brightness diverges for values of (3 < 1/2, so that a cut-off 
radius must be used for the luminosity extrapolation. As neither 
of the two choices of radius used by OP04 is entirely unbiased, 
we performed the Lx /Tx analysis for four choices of luminosity 
cut-off radius. These were: the ROSAT extraction radius, r cut , as 
given in Table EI the fixed over-density radius, rsoo, as calculated 
by OP04; a physical radius of 200 kpc, r p h y s ; and r4 Core , defined as 
4 x the group core radius, r c . We used the luminosity measured to 
r cut and the fitted /3-model parameters to calculate the luminosity 
to rsoo, r p hys and r/t core . We neglected the axial ratio parameter, e, 
included in the OP04 fits, which means that the extrapolated lumi- 
nosities will be slightly overestimated for groups that have a large 
axial ratio. The maximum value of e for a source in our sample is 
2.65 (NGC 4589), but most e values are in the range of 1.0-1.5. 
The radio-loud and radio-quiet groups have similar distributions of 
e, and so our conclusions should not be affected. OP04 could not 
fit a /3 model for five sources, and for these we used (3 = 0.45 (the 
median value measured for the sample) and r c determined by us- 
ing the rough correlation between X-ray luminosity and r c shown 
inFig.Q-Lx = 2.9 x 10 41 (r c /kpc) i n ergs s _1 . 

Another consideration is whether the measured X-ray lumi- 
nosities contain any contribution from sources other than the group 
gas, i.e. from AGN and X-ray binaries. Chandra observations can 
resolve the X-ray binary population and allow the integra ted lu- 
minosity from X-ray binaries to be determined. iKraft et alj i200ll> 
find an integrated X-ray luminosity of ~ 5 x 10 3 9 ergs s -1 for 
the po pulation of X-ray binaries in Centaurus A. IZhang and XvJ 
(2004) examine the X-ray binary population in NGC 1407, one of 
the groups in the sample, which has a comparatively prominent X- 
ray binary population, and find that resolved X-ray binaries account 
for < 20 percent of the galaxy-scale extended X-ray emission de- 
tected with Chandra (or ~ 1.5 x 10 40 ergs s _1 , only ~ 3 per 
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cent of the i?0£Ar-measured Lx for the group-scale emission from 
NGC 1407). We conclude that contamination from X-ray binaries 
is not important to our results. Any contamination from X-ray bi- 
nary emission that cannot be resolved by Chandra should affect the 
radio-quiet and radio-loud groups in the same way. We discuss the 
issue of AGN contamination in Section lTTl 

2.5 L x /T x relations 

For each of the 12 combinations of radio and radius cut-offs, we 
fitted an Lx /Tx relation to the radio-quiet and radio-loud subsam- 
ples. Using the temperature measurements of OP04 and the lumi- 
nosities determined as described above, we fitted the OLS (ordinary 
least squares) bisector fasobe et alll990h to each dataset for consis- 
tency with OP04. Table[4]gives the parameters for the resulting fits. 
We plot Lx vs. Tx for each radio cut-off in Fig.|2|(showing only 
the results for r cu t) with the best-fitting radio-quiet relation over- 
plotted. There is a clear tendency for radio-loud groups to be on 
the hotter side of the radio-quiet relation for all combinations of ra- 
dio cut-off and luminosity extraction radius. (One major exception 
to this trend is NGC 3557, a radio-loud group that is much cooler 
than the prediction for its luminosity, but this may be due to includ- 
ing the cooler galaxy atmospheres of group members, which are 
prominent in a Chandra observation.) 

We tested the significance of the trends illustrated in Fig. [2] 
for each choice of luminosity cut-off. We tranformed the lumi- 
nosity values into a predicted temperature using the appropriate 
best-fitting radio-quiet Lx /Tx relation, as given in Table E| We 
then rotated the coordinate system by —45 degrees so that the x- 
coordinate in temperature corresponds to perpendicular distance 
from the best-fitting line. We then performed a 1-D Kolmogorov- 
Smirnov test comparing the distributions of x (perpendicular dis- 
tance from the line) for the radio-quiet and radio-loud samples of 
each case. The results are given in Tabled] In all but two cases, the 
probability that the two subsamples have the same parent popula- 
tion is < 10 percent, and in more than half of the remaining cases, 
the probability is less than < 5 percent. 

We therefore conclude that there is good evidence that radio- 
loud and radio-quiet groups display different gas properties. The 
choice of radio-luminosity cut-off does not appear to have an im- 
portant effect, whereas the choice of X-ray luminosity radius is cru- 
cial. To obtain a more consistent set of X-ray luminosity and tem- 
perature measurements, higher sensitivity data would be required 
so that the temperature and luminosity could be measured to much 
larger radii, and the need for extrapolation would be reduced. 



3 INTERPRETATION OF THE RESULTS 

The gas properties of radio-loud and radio-quiet groups differ in 
the sense that radio-loud groups of a given luminosity are likely 
to be hotter than the radio-quiet groups. There are several possible 
explanations for this result, which we examine in this section. 

The first question is whether contamination from AGN emis- 
sion could have affected the results for radio-loud groups. We test 
this in Section l3~T1 We then consider three possible physical origins 
for the difference in the properties of radio-loud and radio-quiet 
groups: radio-source heating (Model I), a luminosity deficit caused 
by the radio source (Model II), or an external mechanism that is 
triggering the radio source and heating the gas (Model III). Mod- 
els I and II are related, since any temperature increase would lead 
to an increase in pressure, expansion and a subsequent decrease in 



density on timescales determined by the sound speed; however, as 
a first step in understanding the radio-source impact, it is important 
to test whether the primary effect seen in our results comes from 
an increase in the group temperature or a decrease in X-ray lumi- 
nosity. In Section l3~2l we test Model I (radio-source heating), by 
examining whether there is any evidence that the strength of the in- 
ferred heating of radio-loud groups correlates with the properties of 
the radio sources. In Sections f3.3l and l3~4l we carry out two investi- 
gations to test Model II (a radio-source-induced luminosity decre- 
ment). Firstly, we examine the distribution of gas, parametrised by 
the /3-model, to determine whether this differs for radio-quiet and 
radio-loud groups. We then examine the correlation between Lx 
and optical luminosity to see whether the X-ray luminosities of 
radio-loud groups are lower relative to their optical luminosities 
than is the case for radio-quiet groups, as would be expected in 
the second model. In Section l3~5l we test Model III (an external 
mechanism), by examining the optical properties of the two sub- 
samples to find out whether the radio-loud groups could be in a spe- 
cific evolutionary state, different to that of the radio-quiet groups, 
where the triggering of radio sources and the heating of gas might 
be favoured. Finally, in Section l3~6l we present a study of archive 
Chandra observations of the radio-loud groups to look for further 
evidence of radio-source heating and interactions between the radio 
source and group gas. 

3.1 AGN contamination and the reliability of the OP04 
results 

OP04 state that they have taken into consideration any contribution 
from a central AGN to the X-ray emission via their point-source 
exclusion method. Central AGN were excluded in 22 cases (Os- 
mond, private communication); however, this does not include all 
of the radio-loud groups in our sample. Contamination from strong 
non-thermal emission might result in higher fitted gas temperatures, 
which could shift groups with AGN on the Lx /Tx plane. How- 
ever, if there were a large contribution from AGN-related X-ray 
emission, then the measured X-ray luminosity from the group gas 
would be overestimated; this would act in the opposite sense. It is 
therefore essential to check that AGN contamination is not lead- 
ing to spuriously high temperatures or overestimation of the group 
luminosities, and so we felt it necessary to confirm the results of 
OP04 for several radio-loud groups. In addition, as our results rely 
strongly on the accuracy of OP04's ROSAT analysis, we felt it was 
appropriate to independently measure the temperature and lumi- 
nosity of 3 radio-loud groups, to test for AGN contamination, and 
3 radio-quiet groups, to confirm the reliability of the luminosity and 
temperature measurements. 

From the XMM analysis of ICroston et alJ l2003h . we found 
that the AGN in 3C 66B and 3C 449 (which are more powerful 
than most of the radio sources in the sample we use here) would not 
significantly contaminate the spectral measurements for the group 
atmospheres. The measured temperatures and luminosities in our 
sample are therefore likely to be reliable for groups with Lx > 
10 42 erg s" 1 . For this reason, we selected three groups with lower 
X-ray luminosities: NGC 4261, NGC 1407, where OP04 did not 
exclude the central AGN, and HCG 90. 

In order to test the influence of AGN contamination, we ex- 
tracted a spectrum for the extraction region used by OP04 (a cir- 
cle of radius r cu t), but excluding the central two arcminutes, so 
as to ensure that the majority of the AGN emission was excluded. 
We used a surrounding annulus for background and excluded any 
contaminating point sources by eye. Our analysis is therefore sig- 




Table 4. K-S test results and best-fitting slopes and intercepts for the hx/Tx relations as described in the text. D is the K-S statistic for the given pair of 
radio-quiet and radio-loud samples, and Prob is the null hypothesis probability for obtaining the given value of D. 





Slope 


Intercept 


D 


Prob 


Slope 


Intercept 


D 


Prob 


Slope 


Intercept 


D 


Prob 


r C ut (RL) 


3.30±0.86 


42.15±0.15 


0.488 


0.048 


3.33±0.86 


42.16±0.15 


0.416 


0.114 


2.75±0.71 


42.03±0.18 


0.464 


0.054 


r C ut (RQ) 


3.64±0.61 


42.65±0.14 






3.76±0.72 


42.61 ±0.14 






3.90±0.81 


42.60±0.10 






I"500 (RL) 


3.07±0.72 


42.35±0.12 


0.488 


0.048 


3.19±0.73 


42.35±0.12 


0.434 


0.088 


2.65±0.59 


42.24±0.14 


0.517 


0.022 


rsoo (RQ) 


3.36±0.51 


42.77±0.12 






3.37±0.57 


42.73±0.13 






3.58±0.70 


42.73±0.09 






r phys (RL) 


2.76±0.72 


42.15±0.12 


0.488 


0.048 


2.84±0.72 


42.15±0.13 


0.452 


0.068 


2.22±0.51 


42.01±0.12 


0.536 


0.016 


r ph ys (RQ) 


3.05±0.50 


42.57±0.15 






3.08±0.56 


42.54±0.15 






3.34±0.6B 


42.56±0.10 






r 4c ore (RL) 


3.42±1.10 


41.47±0.1S 


0.450 


0.084 


3.40±1.08 


41.49±0.19 


0.416 


0.114 


2.68±0.91 


41.29±0.20 


0.518 


0.022 


r 4c ore (RQ) 


3.83±0.77 


42.07±0.26 






4.02±0.94 


42.02±0.26 






4.13±1.04 


42.04±0.18 







nificantly cruder than that of OP04, who carried out a more com- 
plicated background estimation and performed good-time interval 
analysis. In all cases our measured temperature is the same within 
the (lcr) errors. However, the measured luminosities in two cases 
are significantly lower than the OP04 results. We conclude that 



there is no risk that the observed difference in the temperature dis- 
tribution of radio-quiet and radio-loud groups is due to spuriously 
high temperatures as a result of AGN contamination. The lower 
luminosity for NGC 4261 is likely to be due to our larger AGN 
exclusion region (which must include significant group emission). 
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In the case of NGC 1407, the slightly lower luminosity is likely to 
be because OP04 did not exclude the AGN. If the AGN contributes 
some of the measured luminosity in a few sources, then the OP04 
luminosity measurements for some of the radio-loud groups may 
be slightly overestimated, which would mean that the significance 
of the effect we observe is underestimated. If the true errors on 
temperature for these groups are slightly larger than those given by 
OP04, this would not have any effect on the K-S test results and our 
conclusions. Finally, since OP04 excluded the AGN in the poorest 
radio-loud groups, including NGC 1052 and IC 1459, we conclude 
that AGN contamination is unlikely to be a problem for our results. 

Even if the AGN emission does not significantly affect the 
measured temperature, it can have a more important effect on the 
surface brightness profile. OP04's neglect of an AGN may affect 
the fitted /3-model parameters for a few groups. Several of the most 
powerful radio galaxies in the sample (e.g. NGC 383, from which 
the AGN was not removed, and NGC 4261) are among the groups 
fitted with a second central /3-model. It is possible that these inner 
/3-models are, in fact, principally modelling the point-source emis- 
sion from the central AGN, although there is also evid ence for a 
galaxy atmosphere in NGC 383 jHardcastle et"ai]l2002h . Since an 
AGN component would be modelled out in this way, the /3-model 
parameters for the extended emission should be reliable. That a sec- 
ond /3-model was not required for many of the radio-loud groups, 
in combination with the lack of any effect on the measured temper- 
atures, supports the conclusion that for most of the less powerful 
radio sources the AGN contribution to the X-ray emission is not 
significant. 

The radio-quiet groups chosen to check the reliability of the 
OP04 results were selected to cover a wide range in X-ray luminos- 
ity and temperature. They are NGC 97, NGC 720, and NGC 4325. 
For each group, we extracted a spectrum using the extraction re- 
gion of OP04 (a circle of radius r cut ) using the same background 
and point-source identification methods as above, and fitted a mekal 
model with free abundance to determine the temperature and X-ray 
flux of the group atmosphere. In all cases the results are in good 
agreement with those of OP04, so we conclude that the OP04 lumi- 
nosity and temperature determinations are reliable both for radio- 
quiet and radio-loud groups. 

3.2 Testing Model I: correlations with radio luminosity 

The results presented in Section|2|strongly suggest that radio galax- 
ies are having an important effect on the properties of the surround- 
ing group gas. We therefore decided to investigate whether there is 
any relationship between the observed "temperature excesses" and 
the radio properties of the associated sources, as would be expected 
in a radio-source heating model. In the following analysis, we use 
only the results for radio-luminosity cut-off cO so as to include the 
widest range of radio powers. 

We used the 1.4-GHz radio luminosity density (Table|2j as a 
measure of the amount of energy a given source might be able to 
provide. Radio luminosity is not an ideal indicator of radio source 
energy input, as the amount of energy transferred from the expand- 
ing radio plasma to the surrounding gas depends on several factors, 
such as source size and age, which are in many cases unknown. 
The size of the source depends on its angle to the line of sight, 
which is usually poorly constrained. The source's age, which is 
also required to estimate the total energy input, is even more dif- 
ficult to determine. Indeed, a few of the radio sources in the sample 
are unresolved with NVSS and FIRST, and do not have identifiable 
double-lobed structure. However, low-frequency radio luminosity 



should trace the jet kinetic luminosity reasonably well, so that in 
the absence of useful information on the sizes and ages of all of 
the sources, the 1.4-GHz luminosity is the best measure of energy 
input available. 

We first compared the radio luminosity with AT, the differ- 
ence between the measured temperature and that predicted by the 
appropriate RQ Lx /Tx relation. Table|5|gives the results of Spear- 
man rank correlation tests for each choice of X-ray luminosity ra- 
dius. For the Spearman test, groups with no temperature increase 
were assigned AT = 0: this applies to 5/19 groups for r cu t and 
)"500, 6/19 for r p h ys , and 3/19 for r core . In all cases there is little 
evidence for a correlation. This is perhaps unsurprising, because 
the observed temperature increase produced by a radio galaxy of 
given luminosity should depend not only on the unknown proper- 
ties of the radio source, as mentioned above, but also on the heat 
capacity of the group gas being heated. A similar analysis using the 
fractional temperature change, AT/T, did not give an improved 
correlation. 

We therefore estimated the heat capacity of each group's en- 
vironment using the spectral and spatial properties of the X-ray 
emission given by OP04. The heat capacity is C — (3/2)Nk, 
where TV is the total number of particles, obtained by integrating 
over the density profile using the best-fitting /3-model parameters, 
with a central proton density obtained from the X-ray luminosity 
and best-fitting mekal model parameters. 

To study the relationship between the observed "heating" and 
radio luminosity, we examined the correlation between E req , the 
energy required to heat the gas in a given group from the pre- 
dicted temperature to the measured temperature (CAT), and L1.4. 
The heat capacities were calculated separately for each of the four 
choices of limiting radius. Fig. [5] shows the relationship between 
L1.4 and E req for each choice of r. For groups with no tempera- 
ture excess, we calculated an upper limit to the energy input by de- 
termining the amount of energy that would be required to shift the 
group significantly to the 'hotter' side of the appropriate Lx /Tx 
relation. As the sample includes upper limits, we used survival 
analysis techniques to determine the generalised Kendall's r corre- 
lation coefficient using ASURV lLavallevetalll992l) . Tablel51con- 
tains the results of the correlation analysis for each case. There is a 
less than 5 percent probability of obtaining the measured value of 
r by chance for 2 out of 4 cases. The high value of 14 per cent for 
r4core is probably because in many cases r4 Core is physically small 
compared to the other choices for r, so that the heat capacity does 
not include much of the gas. For all four choices of radius, there is 
a stronger correlation here than was found for AT alone. 

Since the calculated heat capacity is related to the measured 
X-ray luminosity, we were concerned that the correlation between 
L1.4 and E req could be caused by an Lx / L\a correlation due to 
the flux limits in the X-ray and radio samples. We therefore also 
carried out Spearman rank tests to look for a correlation between 
heat capacity and L1.4. Those results are also included in Table|5] 
and show that the correlation between L1.4 and heat capacity is 
much weaker than that with E req in all cases. 

The presence of a correlation (albeit with a large scatter) be- 
tween radio luminosity and the energy input needed to cause the 
observed temperature increase provides support for a model where 
the temperature increase is due to radio-source heating. The large 
scatter is not surprising, given the many unknown factors, such as 
source size and age, that would affect the amount of observed heat- 
ing. 
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Figure 3. Plot of L1.4, the 1.4-GHz radio luminosity, vs. E re q, the necessary heat in put to produce the temperature increase, for r cu t (top left), rsoo 
(top right), r p ;j ys (bottom left) and r4 core (bottom right). + symbols are groups with temperature excesses; arrows indicate upper limits for groups with no 
observed excess, calculated as described in the text. 



Table 5. Results of correlation analysis for L1.4 and A.T,E req and C. The sample contains 19 groups, so for all of the tests there are 17 degrees of freedom. 



Dataset 




AT 






C 




r a 


Probability 


Kendall's r 


Probability 


r s 


Probability 




0.184 


0.450 


1.992 


0.046 


0.312 


0.193 


^500 


0.214 


0.379 


1.921 


0.055 


0.325 


0.175 


^phys 


0.369 


0.121 


2.555 


0.011 


0.312 


0.193 


T'4core 


0.168 


0.491 


1.472 


0.141 


0.314 


0.190 



3.3 Testing Model II: /3-model properties 

We studied the /3-model properties of the radio-quiet and radio- 
loud groups to determine whether the spatial distribution of gas 
is affected by the presence of a radio source. We compared three 
parameters, Pfu, the fitted value of /3 from OP04, f3 sp ec, the spec- 
troscopic value of f3 from OP04, and Rp = f3 spsc / flfu. For each 
of these properties, we compared the values for the radio-quiet and 
radio-loud subsamples for radio cut-off cO. We performed a 1-D 



K-S test to determine whether the "radio-quiet" and "radio-loud" 
subsamples differed in each case. Fig. |4] shows histograms of the 
distributions of /3/it, f5 apec , and Rp for the radio-quiet and radio- 
loud samples. 

There is no evidence that the parent population differs signif- 
icantly for any of the three parameters. Since the /3 spec have large 
errors, using the K-S test to compare the two samples may not be 
very reliable. We also used a median test to compare the two sam- 
ples, and find a low probability that the distributions of f3fu have a 
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Figure 4. Histograms showing the distribution of fija (top left), f5 spec (top right), and Hp (bottom) for the sample. Filled rectangles are the radio-quiet 
sample, with the radio-loud sample overplotted as hatched rectangles. 



different median. However, there is a probability of ~ 93 percent 
that the distributions of both /3 apec and Hp have different medi- 
ans for radio-quiet and radio-loud groups, in the sense that /3 S p SC 
is higher for RQ groups. This is not a strong result, because of the 
large errors on f3 spe c and therefore Rp. There is no evidence that 
RL groups have flatter profiles than RQ groups, as might be ex- 
pected if the luminosity had been significantly decreased as a result 
of radio-galaxy input. In Section |4] we present further discussion 
of how the group density distribution might be affected by radio- 
galaxy energy input 

3.4 Testing Model II: L x /L B relation 

The X-ray and optical luminosities of groups are correlated, be- 
cause gas mass and galaxy mass should scale similarly. OP04 show 



that such a correlation exists for their sample. If the effects we ob- 
serve in Fig.|2|are caused by a decrease in X-ray luminosity in the 
radio-loud groups, then the Lx/Lb relation should be affected: 
radio-loud groups should have a lower X-ray luminosity relative to 
their optical luminosity. In Fig.|5| we show the Lx I Lb relation for 
radio-quiet and radio-loud groups (using cl and r cut ). Unlike what 
is seen for the Lx /T x relation, there is no apparent difference in 
the two subsamples. We note that that the radio luminosity will be 
related to is, which may introduce a slight bias, but this should 
not affect the X-ray-to-optical luminosity ratio of the radio-loud 
groups. Therefore, this is a strong argument against X-ray lumi- 
nosity decrements in radio-loud groups, as the radio source should 
not affect the optical luminosity of the group. 
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Figure 5. The Lx/Lb relation for radio-quiet and radio-loud groups. 
Symbols are as for previous figures. 



3.5 Testing Model III: optical properties of the RL and RQ 
subsamples 

To test the possibility that radio-quiet and radio-loud groups are in 
different stages of evolution, so that an external mechanism might 
be causing the heating effect, we compared two measures of their 
optical properties. If radio-quiet groups are at a different stage of 
evolution from radio-loud groups, then N ga i s , the number of galax- 
ies in the group, might be expected to differ, in the sense that older 
groups might be expected to have fewer members. Older groups 
might also have a larger ratio between the brightest and second- 
brightest group galaxies, as the largest mergers should have already 
taken place. We therefore compared both N ga is and L12, the lumi- 
nosity ratio of the brightest to second-brightest group galaxy, as de- 
termined by OP04, for the two subsamples (using cl and r cu t). In 
Fig. [6] we show histograms of the distribution of these parameters 
for the two subsamples. In neither case is there a significant differ- 
ence in the distribution for the two subsamples. Both have a peak 
in the value of N ga i a between 5 and 10, and the preferred value 
of I/i2 for both subsamples is < 2, so that most groups have at 
least one reasonably large secondary galaxy, whether radio-loud or 
not. There are no radio-loud groups in the sample with L12 > 10, 
whereas there are two radio-quiet groups with L12 > 20. However, 
a K-S test indicates no significant difference in the distributions 
of either. A median test also finds no significant difference in the 
medians of the two subsamples for either parameter. We conclude 
that the galaxy distributions in radio-quiet and radio-loud elliptical- 
dominated groups are similar, so that there is no evidence that the 
two subsamples are at different evolutionary stages. However, a 
thorough investigation using more sophisticated measures of group 
history is required to test this model fully. 

3.6 Chandra and XMM-Newton observations of heating and 
interactions in the RL groups 

Many of the radio-loud groups in this sample have been observed 
with Chandra or XMM-Newton. The high resolution of Chandra is 
excellent for resolving an AGN component, and for detecting inner 
structure in groups; however, as a result of the high resolution, its 



sensitivity to extended emission is reduced, so that in some cases it 
is unable to detect low surface-brightness emission from the weaker 
groups. In those cases ROSAT temperature and luminosity measure- 
ments are likely to be superior. We discuss here a few groups in the 
sample for which Chandra and XMM-Newton observations show 
evidence for heating or interactions between a radio-source and its 
environment. 

3.6.1 NGC4261 

An XMM-Newton observation of this group was made on 16 De- 
cember 2001 (ObsID 0056340101). An analysis of the extended 
emissi on has not yet been been published, although Sambru na et alJ 
(2003) presented an analysis of the nuclear emission. We ex- 
tracted the archive XMM-Newton data and carried ou t stand ard pro- 
cessing and filtering as described in lCroston et alJ 120031) . Fig. 
shows the adaptively smoothed, background point-source sub- 
tracted, vignetting-corrected 0.5 - 5.0 keV image made from com- 
bined MOS1, MOS2 and pn data, with radio contours from a 1.4- 
GHz map made from VLA archive data overlaid. This figure illus- 
trates a striking relationship between radio and X-ray morphology 
similar to that seen in 3C 66B jCroston et alj|2003l) . It is interest- 
ing that such evidence for interactions between the radio source 
and hot gas on large scales is found in every FRI radio galaxy for 
which deep XMM images of the large-scale structure exist [see also 
lEvansetalJ<2004) l. 

3.6.2 NGC4636 

The Chandra observation of the atmosphere surrounding 
NGC 4636 revealed striking substructure in the form of symmetri- 
cal bright arms Ijones e t al. 2002). Jones et al. find that the leading 
edges of the arms are ~ 30 percent hotter than the surrounding gas, 
and postulate a model in which the arms are produced by shocks 
driven by symmetric off-centre AGN outbursts. The central radio 
source app ears to be extended to the northeast and southwest (e.g. 
Birkinshaw and Davieslll985T) : however, it is small and too weak 
to have produced the shock-heating. Jones et al. argue that this 
indicates that a more direct nuclear outburst may have produced 
the shocks, but it is difficult to think of such a mechanism. 
The radio source could be more extended at low frequencies; 
however, it remains unlikely that a currently active radio source is 
producing the shocks. It is possible that a previous radio outburst 
is responsible, although it is unclear for how long the sharp density 
and temperature structure would persist. The shock-heated arms 
of gas in this group may be the main contributor to the overall 
temperature of 0.84 keV measured by ROSAT. As NGC 4636 has 
one of the largest temperature excesses, this example strongly 
suggests that we are indeed identifying groups with interesting 
AGN/group interactions. 

3.6.3 NGC 1052 

The Chandra observation is too short to detect much low surface- 
brightness emission from this poor group. However, as shown 
in Fig. [8] th ere is evidence for radio-related X-ray emission, as 
discussed by Kadler et al. (2004). They attribute most of the X- 
ray emission to the jet; however, the distribution of X-ray counts 
around the eastern radio lobe seems to be reminiscent of the bright 
shell of hot gas around the southwestern inner lobe of Centaurus A. 
The two systems are remarkably similar: they both consist of small, 
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Figure 7. An adaptively smoothed, background point-source subtracted, vignetting-corrected 0.5 - 5.0 keV image of NGC 4261 made from the combined 
MOS 1 , MOS2 and pn data from the archive XMM observation described in the text. Clear evidence for interactions between the radio galaxy (3C 270) and gas 
environment are seen in the form of holes in the X-ray surface brightness at the positions of both radio lobes. 



double-lobed radio sources, which are likely to be young or restart- 
ing. They have galaxy atmospheres of similar X-ray luminosities 
and temperatures. These Chandra observations therefore suggest 
that the young radio source in NGC 1052 may be shock-heating 
its surroundings. The measured X-ray temperature from ROSAT 
may contain a large contribution from these radio-related regions, 
although it is also possible that the entire environment has been 



heated, as we have argued to be the case for 3C 66B ICroston et alj 
2003). 

3.6.4 HCG62 

The X-ray emission from HCG 62 provides one of the c learest 
examples of "holes" in a group atmosphere <VrtileketalJl200d) . 
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Figure 8. Gaussian smoothed image of the Chandra data for NGC 1052 
in the energy range 0.5 - 5.0 keV. Radio contours are overlaid from a map 
made from archive VLA data IWrobel 19841) . 



However, the current AGN is a weak radio emitter, and does not 
show any extension. As with NGC 4636, it is plausible that a previ- 
ous radio outburst (which may still be detectable in low-frequency 
radio observations) has produced the observed X-ray structure. 

3.6.5 NGC 5044 

This group shows prominent substructure in the Chandra image of 
iBuote etail 1200 3h . who associate holes in the gas with the radio 
source. Although the NVSS image of this source does not show 
evidence for any extended radio emission, they suggest that obser- 
vations at lower frequency might reveal the presence of radio emis- 
sion filling the cavities. We have been unable to resolve the AGN or 
detect larger-scale radio emission in our analysis of VLA archive 
data. 

3. 6. 6 Summary of the Chandra and XMM-Newton observations 

The Chandra and XMM-Newton observations discussed above 
show that radio-source/group interactions are complex. In two 
cases, NGC 1052 and NGC 4261, there is clear evidence that the X- 
ray structure has been affected by the current radio galaxies. In two 
groups where there are no detected large-scale radio lobes, HCG 62 
and NGC 4636, the Chandra observations reveal striking X-ray 
morphologies suggestive of outbursts from the AGN. Smaller-scale 
substructure is also present in NGC 5044. Finally, localised heating 
appears to be present in NGC 1052 and NGC 4636, suggesting that 
the heating effects we observe from the ROSAT sample could be 
caused by several different processes. 



4 EVIDENCE FOR RADIO-SOURCE HEATING? 

The results presented in Section |2| support the argument that ra- 
dio galaxies have an important effect on the X-ray properties, and 
therefore the physi cal conditions, of g roup gas, as suggested by 
previous work ("e.g. lCroston et all2003h . This is shown by the dif- 
ference in Lx /Tx properties of the two subsamples, and by the 
high incidence of radio-related substructure in radio-loud groups. 

In the next two sections we discuss the three models of Sec- 
tion|5|in detail. In Section l4~T1 we compare models I and II, and ar- 
gue in favour of a radio-source heating model (Model I) and against 
a luminosity deficit (Model II). In Section l4~2l we consider one 
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possible external mechanism that could lead to a Model III expla- 
nation, that of mergers and interactions, and argue against such a 
model. In Section l4~3l we attempt to explain the results for all of the 
radio-loud groups in the context of a model of radio-source heating 
and discuss what can be inferred about the heating processes. 

4.1 Temperature increase vs. luminosity decrement 

Radio galaxies must displace large amounts of gas and this could 
have a significant effect on their luminosity. For 3C 66B, which 
is larger than most of the sources in the samples studied here, we 
calculated that the gas with which the radio source can have di- 
rectly interacted provides only 7 percent of the group's luminosity 
ICrostonet afl Eo03). It is therefore unlikely that removal of gas 
by the radio galaxy could produce the luminosity deficits needed 
by this model, in some cases an order of magnitude in luminosity. 
However, the group luminosity will also be decreased if a signif- 
icant fraction of the jet kinetic energy is transferred into potential 
energy. 

In the context of preheating mod els of energy injection into 
group gas, it has been ar gued (e.g. iMetzler and Evrarcj 1 19941 
Helsdon and Ponman 2000) that the main effect of the energy in- 
jection will be an increase in the group's potential energy, so that 
the central density decreases (and hence luminosity will decrease). 
While (by the virial theorem) this will eventually be the case, heat- 
ing effects are still l ikely to be detectable on shorter timescales. 
More recently, iKavl 120041) has carried out cosmological simula- 
tions of cluster formation including cooling and feedback (which 
could be due to AGN or a different energy source such as supernova 
winds) and find gas properties in agreement with observations, with 
a ~ 10 percent increase in temperature at the virial radius. Their 
simulations consider only massive clusters, but suggest that at least 
some fraction of AGN-injected energy is likely to end up in the 
thermal energy of the group. 

In some of the groups in our sample an order of magnitude 
decrease in luminosity is required: the density would have to be 
dramatically reduced to produce such an effect. A strong argument 
against such large luminosity deficits in radio-loud groups comes 
from the Lx / L b relation discussed in Section l3~4l We find that the 
radio-loud groups follow the same trend as the radio-quiet groups 
and show no evidence for having lower X-ray luminosities relative 
to their optical luminosities, as would be expected if an X-ray lu- 
minosity decrease had been caused by the radio galaxy. 

Another strong argument in favour of heating as the dominant 
effect, as opposed to a change in luminosity, is the result of Sec- 
tion l3.2l where we found evidence for a correlation between radio 
luminosity and the energy required to heat the gas from the pre- 
dicted to the observed temperature. This result would be harder to 
explain in a model where the radio-source's impact was principally 
on the group's luminosity. 

As shown in Section l3~6l several groups with a temperature ex- 
cess possess additional evidence for radio-source heating. In Fig.|9| 
we show the Lx /Tx relation for cl and r cut , with NGC 4636 
and NGC 1052, as well as 3C 66B 3C 449 and NGC 6251 
ICroston et afll2003b lEvans et alJl2004h marked to illustrate how 
they compare to the sample studied here. We conclude, based on 
the additional evidence for heating in several sources, and the ar- 
guments above, that heating is a more plausible explanation than a 
radio-source induced luminosity deficit. 

On the longest timescales, the energy injected by radio galax- 
ies into group or cluster gas must predominantly end up as poten- 
tial energy, and any long-term temperature increase will be small. 
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Figure 9. Lx/Tx relation (cl, r cu t) showing the po sitions of the 
three X M M-observed radio-gal axy environments presented in Cros ton et alj 
i 2003) and Evans et al. 1 2004), as well as the two OP04 groups with addi- 
tional strong evidence for heating, NGC 1052 and NGC 4636, marked. + 
symbols indicate radio-quiet groups, and hollow squares radio-loud groups. 



However, information about the energy injection cannot travel 
faster than the sound speed, so that temperature effects may be 
detectable for a few sound-crossing times. That an Lx/Tx rela- 
tion for elliptical-dominated groups exists at all is evidence that, 
on average, the temperature increase must disappear on timescales 
less than the radio-source recurrence time; for a 50 per cent radio- 
galaxy duty cycle, this is comparable to a few sound crossing times 
in a typical group. It is plausible that occasionally a second AGN 
outburst would occur before the group has recovered from the pre- 
vious outburst, so that the temperature increase is more persistent; 
3C 66B could be one example of such a system. Our results there- 
fore suggest that we are detecting the short-term effects of radio- 
source heating in many elliptical-dominated groups. The fact that 
we have found no systematic differences in the properties of radio- 
loud and radio-quiet groups, other than their Lx /Tx relations, is 
consistent with a model in which all elliptical-dominated groups 
have had similar numbers of radio-galaxy outbursts averaged over 
the group lifetime, affecting the groups by causing a temporary in- 
crease in temperature, with less easily detectable long-term effects 
on the group luminosity and surface brightness distribution. 

4.2 The evolution of groups: a common cause for heating and 
radio activity? 

An alternative explanation is that some common cause triggers a 
radio source and heats the gas. It is possible that the elliptical- 
dominated groups containing radio sources exist at a particular 
stage in the evolutionary process for groups, where the gas is hotter 
relative to the group luminosity. One possibility is that the inci- 
dence of mergers, and/or the type of mergers that the two subsam- 
ples of groups have undergone, is different for the two subsamples. 

Recent simulations by Rowlev et al. (2004) find that in ma- 
jor mergers (defined as increasing the cluster mass by > 20 per- 
cent) clusters become brighter and heat up roughly parallel to the 
Lx/Tx relation, whereas in minor mergers, the temperature in- 
creases and the luminosity decreases. Although these simulations 



are for more massive clusters, this suggests a possible interpreta- 
tion of our results. If the radio-loud groups had recently undergone, 
or are continuing to undergo, minor mergers, then their Lx /Tx 
properties could be explained. In this model, the radio-quiet groups 
must either have undergone mainly major mergers, or be in more 
isolated environments where mergers are less common. One possi- 
bility is that radio-quiet groups could be more evolved, so that most 
merging has already occurred, and the gas has cooled back to the 
Lx /Tx line. In this scenario, it might be expected that the radio- 
quiet groups would resemble fossil groups where mergers are no 
longer common. 

Such a model of group evolution should be testable using the 
optical properties of the group. However, we found no evidence 
for a difference in either the number of galaxies in the group or 
the degree of dominance by the central galaxy for the radio-quiet 
and radio-loud groups (Section |3.5> . A model of this sort cannot 
be ruled out, as mergers remain the most plausible model for how 
radio galaxies are triggered. However, there does not appear to be 
any evidence that the radio-quiet and radio-loud groups are in dif- 
ferent stages of evolution, based on the comparison of the optical 
group properties. We therefore conclude that radio-source heating 
is the most plausible of the three explanations for our results. 

4.3 Models for radio-source heating 

If the radio-source heating model is adopted, then it is necessary to 
consider whether it is possible to explain the observed results for all 
of the radio-loud groups in the context of this model. In Fig. [2] it is 
apparent that there are some "radio-loud" groups that do not show 
a temperature excess, and one particularly anomalous group that is 
much cooler than predicted (NGC 3557). In addition to the sources 
with no large temperature excess, most of the currently active ra- 
dio sources in the sample are unlikely to be capable of producing 
all of the heating that is observed. The two groups HCG 62 and 
NGC 4636, which show strong evidence for radio-related struc- 
ture and a large heating effect (especially in the latter), are par- 
ticularly problematic. Either there is low-frequency radio structure 
that has not yet been observed, which indicates active or recently 
switched-off large-scale radio jets and lobes, or else the heating ef- 
fects are long-lived. Some of the most powerful radio sources, such 
as 3C 66B, 3C 449, and NGC 4261, are probably capabl e of pro- 
ducing the heating that is observed ICroston et alJl2003h . but the 
heating in many of the radio-loud groups must be longer-lived than 
the radio source. 

It is also interesting to consider the mec hanisms for he ating in 
different stages of radio-source evolution. In lCroston et alj l2003h . 
we argued that large, powerful FRIs are subsonic and likely to be 
heating their surroundings gently via PdV work as the lobes ex- 
pand. However, in the early stages of FRI evolution, the sources 
are known to be overpressured (even assuming minimum energy 
pressures), and the recent observation of a heated shell around the 
inner lobe of Cen A shows that shock-heating is not only likely to 
be an important mechanism in FRIIs, but plays a role in the early 
stages of FRI evolution as well. This process may also be occurring 
in one of the groups in this sample, NGC 1052 (see Section l3.6.3> . 
Finally, the Chandra observations of NGC 4636 show that addi- 
tional mechanisms for AGN heating may exist, as it is difficult to 
explain the morphology of the shocked arms of gas via radio-lobe 
expansion. We conclude that the heating effects found in the study 
of the ROSAT sample presented here are likely to be the result of 
different types of radio-source heating, so that one simple model for 
the entire sample is unlikely to be correct. In some sources there 
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may be small regions of shocked gas, unresolved in the ROSAT 
data, leading to the temperature increases, whereas in others more 
widespread heating is necessary. A detailed analysis of XMM ob- 
servations, which now exist for a significant number of the radio- 
loud groups, would help to investigate these possibilities, as would 
a low-frequency radio study to constrain the properties of the radio 
sources. 

4.4 The importance of radio galaxies in structure formation 
models 

We have shown that radio galaxies at some stage of development 
are present in up to 50 percent of elliptical-dominated groups (Sec- 
tion l2.3t . and have also presented evidence that radio-source heat- 
ing is common. It is therefore of interest to consider whether their 
energy input is of significance in the context of structure forma- 
tion models. We carried out some simple calculations to deter- 
mine whether the energy input from low-power radio sources in 
elliptical-dominated groups could be important in the context of 
the energy-injection requirements of structure formation models. 

We estimated the energy input rate from the average radio 
source in the sample by taking the average value of L1.4 for the 
19 radio-loud groups (3.5 x 10 23 W H z -1 ) and scaling the kin etic 
luminosity of 3C 3 1 (from the model o flLaineand Bridlel2002h bv 
the ratio of L1.4 for the average radio source and 3C 31, which 
gives 7 x 10 35 W. We then assumed that 1/3 of the kinetic lumi- 
nosity is transferred to the group gas, a conservative lower limit. 
We assumed a duty cycle of 50 percent, based on the fraction of 
radio sources in the Garcia catalogue (Section |2.3> . This gives a 
typical energy input rate of 2.3 x 10 51 keV/s over the lifetime of 
the group. The energy rate per particle was determined using the 
average number of particles in the group (determined as part of 
the heat capacity calculations in Section l3~2l to be 6.6 x 10~ 18 
keV/particle/s. Assuming the injection energy is required to be of 
order 1 keV/particle (e.g. IWu et al]l200 Q1. then the average radio 
source in the sample can provide the necessary energy input over 
~ 5 x 10 9 years, which is a plausible group lifetime (this is about 
10 times the standard radio-galaxy lifetime, so that the radio source 
must have ~ 5 active phases during the group lifetime). We con- 
clude that low-power radio sources may be capable of providing the 
necessary energy input in elliptical-dominated groups. As it is the 
elliptical-dominated groups that principally determine the Lx /Tx 
relation for groups, since they dominate the population of groups 
with luminous, group-scale X-ray environments, it is therefore pos- 
sible that the energy input from low-power radio galaxies can ex- 
plain the X-ray properties of groups. 

Other workers have carried out calculations of the energy in- 
put of radio galaxies into clusters. It seems likely that low-power 
radio galaxies can provide sufficient energy on average to balance 
cooling flows (e.g. lFabian et alj2"00 3V however, to explain the clus- 
ter Lx /Tx relati on may req uire more energy than can be supplied 
by FRIs [although Rovchowdhu rv et aljj2004j) found that efferves- 
cent heating by rising bubbles could solve the entropy problem by 
heating at large radii]. The energy contribution from FRII radio 
galaxies and quasars is also likely to be important in the rarer situ- 
ations where they occur. 



5 CONCLUSIONS 

We have presented a detailed study of the gas properties of radio- 
quiet and radio-loud elliptical-dominated groups based on the 



ROSAT groups sample of Osmond & Ponman (2004). We reach 
the following conclusions: 

• 63 percent of the elliptical-dominated groups (19/30) in the 
OP04 sample have associated radio sources at the centre of a dom- 
inant group galaxy. 

• Our sample of elliptical-dominated groups is not significantly 
biased in its radio-loud fraction: the true fraction in the parent pop- 
ulation may be ~ 40 — 50 per cent. 

• Radio-loud groups are likely to have a higher temperature than 
radio-quiet groups of the same luminosity. 

• The energy required to produce the observed temperature ex- 
cess correlates weakly with the 1.4-GHz radio luminosity of the 
sources. 

• The difference in gas properties for radio-quiet and radio-loud 
groups is most plausibly interpreted as evidence for radio-source 
heating. 

• Evidence for radio-source interactions with the surrounding 
gas is found in Chandra or XMM-Newton observations of many of 
the radio-loud groups, although there are also several groups that 
show disturbances not directly related to observable radio structure. 

• The radio-loud groups are at different stages in the heating 
process, so that some may be experiencing shock-heating by young 
radio sources, some are being gently heated by a currently active 
large radio galaxy, and some show longer-lived heating effects from 
a previous generation of radio-source activity. 
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